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ABSTRACT 

We determine the Minkowski functionals for a sample of Abell/ACO clusters, 
401 with measured and 16 with estimated redshifts. The four Minkowski func- 
tionals (including the void probability function and the mean genus) deliver 
a global description of the spatial distribution of clusters on scales from 10 
to 60/i Mpc with a clear geometric interpretation. Comparisons with mock 
catalogues of N-body simulations using different variants of the CDM model 
demonstrate the discriminative power of the description. The standard CDM 
model and the model with tilted perturbation spectrum cannot generate the 
Minkowski functionals of the cluster data, while a model with a cosmological 
constant and a model with breaking of the scale invariance of perturbations 
(BSI) yield compatible results. 

Key words: cosmology - large scale structure - statistical methods - struc- 
ture formation 



1 INTRODUCTION 

The distribution of galaxy clusters has been used for 
a long time as a useful tracer of the large-scale struc- 
ture of the Universe and as a constraint for cosmological 
models. To present knowledge, clusters of galaxies rep- 
resent the largest gravitationally bound entities in the 
hierarchy of cosmic structures. In addition, their ampli- 
fied clustering, with respect to that of galaxies, can be 
reliably measured on large scales (> 20 h~ Mpc) where 
the gravitational evolution still keeps track of the initial 
fluctuation spectrum. 

The first quantitative analyses of the cluster distri- 
bution date back to the pioneering works by Bahcall & 
Soneira (1983) and Klypin & Kopylov (1983) and were 
based on the estimate of the cluster 2-point correlation 
function, £ C c(r). Although dealing with a rather limited 
number of objects, these analyses clearly showed that 
Abell (1958) clusters have a clustering scale length ro, 
defined by £ cc (ro) = 1, which is much larger than that 
of galaxies. Considerable observational effort has been 
made to improve the large-scale mapping of the clus- 
ter distribution. Analyses based on the extension of the 
Abell sample to the southern hemisphere (Abell, Cor- 
win & Olowin 1989, hereafter ACO) and on progres- 



sively larger redshift compilations of Abell and ACO 
clusters (e.g. Postman, Geller & Huchra 1992) basically 
confirmed the original result that ro — 20ft~ 1 Mpc (see 
also Cappi & Maurogordato 1992; Plionis, Valdarnini & 
Jing 1992; Peacock & West 1992). On the other hand, 
new samples based on more objective cluster identifi- 
cation criteria, both in the optical (Dalton et al. 1994; 
Collins et al. 1994) and in the X-ray (e.g. Nichol, Briel 
& Henry 1994; Romer et al. 1994) bands, revealed a 
smaller length scale, ro — 13. . . 16/i~ 1 Mpc. This differ- 
ence might be ascribed to richness contamination from 
projection effects of the Abell/ACO samples (Suther- 
land 1988, see also Jing, Plionis & Valdarnini 1992) or 
to the richness-clustering dependence (Bahcall & West 
1992). In any case, the standard CDM model fails to 
predict the large range of positive cluster correlation 
function (£(r) > up to scales of 70fo -1 Mpc, see White 
et al. 1987, Olivier et al. 1993). 

Although several further analyses have been de- 
voted to compare data on the cluster distribution to the 
predictions of a variety of DM models (e.g. Bahcall & 
Cen 1992; Dalton et al. 1994; Borgani et al. 1995), most 
of them have been mainly based on low-order statistics. 
A more complete description of the cluster distribution 
requires the use of higher-order statistics, which are ca- 
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pable of capturing global aspects of the clustering. One 
example is the genus statistics applied to the cluster 
data by Rhoads, Gott & Postman (1994) and by Pear- 
son et al. (1996). 

A direct measurement of properties of a general 
point distribution can be made using Minkowski func- 
tionals (Mecke, Buchert & Wagner 1994). These mea- 
sures provide a morphological (geometrical and topo- 
logical) description of the point distribution and com- 
prise information about correlation functions of arbi- 
trary order. Application to clusters of galaxies is espe- 
cially promising since clusters bear direct contact to the 
formation of the largest structures in the universe. As 
geometrical characteristics of structure, the Minkowski 
functionals combine both the advantage of intuitive in- 
terpretation of their meaning and the advantage of de- 
livering a quantitative measure. The aim of this paper 
is to apply the Minkowski functional statistics in order 
to compare an extended redshift sample of Abell/ACO 
clusters with mock catalogues obtained from large PM 
N-body simulations, starting from different models for 
the initial power spectrum. To this purpose, we extract 
artificial cluster samples from the simulations, thereby 
reproducing the main characteristics of the observa- 
tional data set (cluster number density, boundary ge- 
ometry, selection functions for galactic absorption and 
luminosity effects). 

Simulations are run for the Standard Cold Dark 
Matter (SCDM) model, a 'Tilted' CDM model (TCDM) 
with v = 0.9 for the post-inflationary spectral index, a 
low-density CDM model with fio = 0.35 and a cos- 
mological constant (ACDM), and a double inflation 
model with Broken Scale Invariance (BSI) of the pri- 
mordial perturbation spectrum. It is generally accepted 
that the SCDM spectrum fails to reproduce the clus- 
ter correlation amplitude and, once normalized on large 
scales to match the CMB anisotropies, it largely over- 
produces clusters. These shortcomings are interpreted 
as the consequence of the wrong shape of the SCDM 
spectrum (see, e.g., Borgani et al. 1996; Mo, Jing & 
White 1996). As a possible remedy, ACDM modifica- 
tions of the SCDM model suitably change the spectrum 
shape at the cost of introducing one additional param- 
eter in comparison to SCDM, while the TCDM model 
and the BSI model change the shape of the power spec- 
trum by introducing one and two additional parameters, 
respectively. 

A brief description of the observational data is pro- 
vided in Section 2. In Section 3 we present the cosmolog- 
ical models on which simulations are based and sketch 
the procedure for cluster selection and mock catalogue 
construction. We define and discuss in Section 4 the 
properties of Minkowski functionals, relations to other 
statistics, the method of analysis and its results. In Sec- 
tion 5 we summarize and draw our conclusions. 



2 OBSERVATIONAL DATA 

We consider an extended redshift sample of Abell and 
ACO clusters with richness R > (Abell 1958; Abell, 
Corwin & Olowin 1989). In the following we will pro- 



vide only a brief description of this sample since more 
details are given in Borgani et al. (1996; see also Plionis 
& Valdarnini 1995). The northern (Abell) part of the 
sample, with declination 8 > —17°, is defined by those 
clusters that have measured redshift z < 0.1, while the 
southern ACO part, with 8 < —17°, is defined by those 
clusters with mio < 17, where mio is the magnitude of 
the tenth brightest cluster galaxy in the magnitude sys- 
tem corrected according to Plionis & Valdarnini (1991). 

The effect of galactic absorption is modelled accord- 
ing to the standard cosecant dependence on the galactic 
latitude b, 

P(|b|) = io Q(1 - csc|b|) , (1) 

with a « 0.3 for the Abell sample (Bahcall & Soneira 
1983; Postman et al. 1989) and a w 0.2 for the ACO 
sample (Batuski et al. 1989). In order to limit the ef- 
fects of galactic absorption we only use clusters with 
|6| > 30°. 

The cluster-redshift selection function P(z) is de- 
termined by fitting the cluster density as a function of z: 

n*) = { \ , i y ! 2 - Zc ' (2) 

w \ A exp(-z/z a ) if z > z c , w 

where A — exp(z c /z ), and z c is the redshift be- 
low which the spatial density of clusters remains con- 
stant (volume-limited sample). We find z c ~ 0.078, 
z « 0.012 and z c « 0.068, z a « 0.014 for Abell and 
ACO samples, respectively. Since the exponential de- 
crease of P(z) can introduce considerable shot noise er- 
rors at large redshifts, we prefer to limit our analysis 
to r max = 240ft _1 Mpc, where the catalogue is approxi- 
mately volume limited. 

There are in total 417 Abell/ACO clusters fulfill- 
ing the above criteria: 262 Abell clusters with mea- 
sured redshifts, and 155 ACO clusters, 16 of which 
have z estimated from the m\o-z relation calibrated 
by Plionis & Valdarnini (1991). These numbers corre- 
spond to (n) A beii = (1.6 ± 0.25) x 10 -5 (h^Mpc)" 3 
and (n) AC o = (2.3 ± 0.3) x 10~ 5 (/i^Mpc)" 3 , for 
the Abell and ACO cluster number densities, respec- 
tively, once corrected for galactic absorbtion according 
to eq. (Q). The density difference is mostly spurious, due 
to the higher sensitivity of the Illa-J emulsion plates 
on which the ACO survey is based. The above den- 
sity values correspond to average cluster separations of 
(d A beii) ^ 40 /i _1 Mpc and (d A co) ^ 35 /i _1 Mpc. 



3 MOCK SAMPLES FROM N-BODY 
SIMULATIONS 

3.1 Cosmological models 

We simulated the evolution of large-scale structure in 
four different types of spatially flat Cold Dark Matter 
models. They are listed below. 

(a) The Standard CDM (SCDM) model with fi = 1 
and h = 0.5 for the Hubble constant in units of 
100 kms -1 Mpc" 1 , which we take as a standard of ref- 
erence. 
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(b) A tilted CDM (TCDM) model, with v = 0.9 for 
the post-inflationary spectral index as arising, for in- 
stance, from power law inflation. 

(c) A low-density CDM (ACDM) model, with 
fio = 0.35, h — 0.7 and spatial flatness restored by a 
cosmological constant term = 0.65. 

(d) A Broken Scale-Invariance (BSI) CDM model, 
which exhibits a step-like primordial spectrum due to 
a double-inflationary scenario. The modification of the 
power spectrum is specified by two parameters, the step 
location at fc^reak = 1.5 h~ Mpc and its relative height 
A = 3. With this parameter choice, the model provides 
a good fit to a variety of observational data (Gottlober, 
Miicket & Starobinsky 1994, Amendola et al. 1995, 
Kates et al. 1995, Ghigna et al. 1996). 

For the CDM model we use the transfer function as 
parameterized by Bardeen et al. (1986), which assumes 
a vanishing baryon contribution. Although the presence 
of baryons is relevant as far as small-scale properties 
(< 5/i _1 Mpc) are concerned, it has negligible influence 
on the large-scale clustering 10/i _1 Mpc), which we 
want to investigate using galaxy clusters. We summarize 
the main model parameters in Table |l| The power spec- 
tra are normalized according to the two year COBE- 
DMR 53 and 90 GHz galactic sky maps following the 
prescription of Gorski et al. (1994) for the Q — 1 models 
and of Stompor, Gorski & Banday (1995) for the ACDM 
model. 

We did not take into account any possible gravita- 
tional wave contributions for the BSI and TCDM model. 
We have also done our analysis for the one year COBE 
normalization and verified that the statistics of the clus- 
ter distribution are insensitive to this reduction of the 
spectrum amplitude (see also Croft & Efstathiou 1994; 
Borgani et al. 1995). Therefore, any uncertainty either 
in our normalization procedure or in the measured level 
of CMB temperature anisotropy has no effect on the fi- 
nal results of the analysis. In particular, our results will 
not change when adopting a power spectrum normaliza- 
tion compatible with the 4 year COBE data (see, e.g., 
Bennett et al. 1996). This normalization would imply 
a reduction of the fluctuation amplitude <rg by about 
11 % for all the models under consideration. 

The spectra we used in our simulations are shown 
in Figure j], where the horizontal bar selects the part 
of the spectrum between the inverse box size and the 
Nyquist frequency which is realized in the simulation. 

We evolve the initial density field starting from 
redshift z — 25 until the present epoch employing a 
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body scheme with 

12; 



300* parti- 



standard PM 

cles (m pa r = 1.3 x 1O 12 /i _1 M ) and N g = 600 3 grid 
cells in a simulation box of L — 500/i _1 Mpc comov- 
ing length a side. This provides a spatial resolution of 
~ 1.7/i _1 Mpc (two cells). We suppose that the simula- 
tion box is large enough to contain all fluctuation modes 
which contribute to the large-scale cluster clumping. In 
order to account for the effect of statistical variance, 
we carried out simulations for four random realizations 
for SCDM and three for ACDM, respectively, while one 
realization was done for TCDM and BSI. 




0.010 0/00 
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Figure 1. The SCDM, TCDM, ACDM and BSI power spec- 
tra. The bar selects the wavenumber range corresponding to 
the adopted box size of 500?t _1 Mpc. 



Table 1. Summary of model parameters 



Model 


h 




Q A 


V 


<T8 


SCDM 


0.5 


1 





1 


1.37 


TCDM 


0.5 


1 





0.9 


1.25 


ACDM 


0.7 


0.35 


0.65 


1 


1.30 


BSI 


0.5 


1 







0.60 



3.2 Mock cluster samples 

We identify clusters from the final particle distribution 
by applying an iterative procedure to the high peaks of 
the density field reconstructed on the 600 3 mesh points 
(see also Klypin & Rhee 1994) . After assigning the den- 
sity field on the grid through a cloud-in-cell interpolat- 
ing scheme, we identify those points which correspond to 
local density maxima. Afterwards, we center a sphere of 
radius 1.5/i _1 Mpc, corresponding to the Abell radius, 
on each of these local maxima and compute the cen- 
ter of mass position for all the particles falling within 
that sphere. This position is used as the new cluster 
center and the procedure is repeated until convergence. 
We find that in general only a few (about 5) iterations 
are required for the cluster coordinates and masses to 
converge to their final values. 

From the resulting list of candidate clusters, we se- 
lect the Nd most massive objects and identify them 
with Abell/ACO clusters. By definition, N ci = (L/d ci ) 3 
is the expected number of clusters within the simulation 
box, having mean separation d c i. 

In order to check the robustness of our cluster 
identification scheme, we also based their identification 
on the friend-of-friend algorithm, with linking length 
b — 0.2 times the mean particle separation (e.g. Frenk et 
al. 1988) . This value of b defines groups bounded by an 
isodensity surface of about 125 times the average den- 
sity and, therefore, to an average internal overdensity 
of about 180 if an isothermal density profile is assumed 
(e.g. Lacey & Cole 1994). Such an overdensity is very 
close to the value expected for a virialized structure re- 
sulting from the spherical top-hat collapse (e.g. Peebles 
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Figure 2. Projected positions of three mock samples in the 
simulation box, showing their overlap. 



Table 2. The Minkowski functionals in three— dimensional 
space expressed in terms of the corresponding geometric 
quantities. 



geometric quantity fi M„ 



V volume V 

A surface 1 A/8 

H integral mean curvature 2 H/2ir 2 

X Euler characteristic 3 3x/4-7r 



1980). We will show elsewhere that the cluster distri- 
butions obtained with these two procedures are very 
similar not only in a statistical sense, but also as far as 
the point-to-point comparison of the cluster positions 
is concerned. 

After generating the cluster distribution in the 
simulation box, we extract mock samples which re- 
produce the same observational features as the real 
Abell/ACO sample. In each box we locate 8 observers 
along the main diagonal axes, each having a distance of 
L/4 = 125/i -1 Mpc from the three closest faces. First 
we include all clusters up to a maximum distance of 
240/i -1 Mpc in each mock sample (see Figure |^, the 
simulation box has periodic boundaries). We then ran- 
domly sample them to get a density distribution repro- 
ducing the observational selection functions for galac- 
tic absorption and redshift extinction (see also Borgani 
et al. 1996). In order to minimize the overlap between 
mock samples, the coordinate systems for two adjacent 
observers are chosen so that the corresponding galac- 
tic planes are orthogonal to each other. Even with this 
choice, it turns out that different mock samples involve 
overlapping volumes and, therefore, they cannot be con- 
sidered as completely independent. 

Minkowski functionals are sensitive to the number 
density. In the case of a Poisson process one may eas- 
ily derive scaling relations with the number density and 
radius using the analytical results of Mecke & Wagner 
(1991) and the homogeneity property of the Minkowski 
functionals, M M (AA r ) = X d ~"M„(A r ), (fi = 0, . . . , 3) 
with a positive real scaling factor A (see Section ^ for 
the definition of the functionals M M of A r produced by 
balls B r of radius r around each cluster). For generic 
distributions no scaling relations are available since the 
scaling properties depend, even without correlations, on 
the dimensionality of the support of the point process. 
Therefore, care has to be taken to reproduce in the sim- 
ulated samples the correct number of Abell and ACO 
clusters separately. To this purpose, after generating a 
mock sample, we randomly degrade the cluster num- 
ber density in the Abell part until the number density 
reaches a fraction (n)Abcii/(w)ACO of the number den- 
sity of the ACO part. Since the overall number density is 
fixed in the whole simulation box, different samples may 
contain different numbers of clusters, with fluctuations 
around 10% and deviations of single samples as large as 
30%. Instead of forcing all the samples to have the same 
number of clusters as the real one, we prefer to main- 




M^A) = M„(gA) 




Mfj_(Bi U B 2 ) = M„(B0 + (£?2 ) - M M (Bj n B 2 ) 




M^d) — > M M (C) as d — -» C 



Figure 3. Minkowski functionals are unique descriptors 
of morphology under the requirements of Motion Invariance 
(A), Additivity (B) and Continuity (C). 

tain such fluctuations and consider them as an effect of 
cosmic variance. We only take care that, after averag- 
ing over the eight observers in each box, the resulting 
average cluster number per sample reproduces the ob- 
servational one. We find that choosing a mean cluster 
separation d c i = 36 /i _1 Mpc, which roughly corresponds 
to the average separation of ACO clusters, always pro- 
duces an average number of Abell and ACO clusters per 
sample which differs by ^ 5% from the corresponding 
numbers of the observational case. 

As we will see, uncertainties in the analysis of 
the Minkowski functionals are largely dominated by 
the observer-to-observer scatter within the same box, 
which is due to the fluctuations in the cluster number 
density. The box-to-box variance of different simula- 
tions, which accounts for the sampling of independent 
patches of the Universe is significantly smaller. 



4 ANALYSIS USING MINKOWSKI 
FUNCTIONALS 

4.1 Definition and general properties 

Let us consider the set of points supplied by cluster 
positions in three-dimensional space. We decorate each 
point with a ball of radius r, thereby creating connec- 
tions between neighbouring balls. We now wish to in- 
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vestigate how the global morphology of the union set of 
these balls changes with the radius r, which is employed 
as a (single) diagnostic parameter. To achieve this, we 
need quantitative measures of geometry and topology 
for bodies in three-dimensional space. 

It seems sensible to request that such measures be 
motion invariant valuations of the bodies, i.e. scalar 
functionals satisfying additivity and invariance under 
rotations and translations, as well as a continuity re- 
quirement (see Figure ^[). The theorem of Hadwiger 
(Hadwiger 1957) tells us that in three dimensions any 
functional satisfying these requirements is a linear com- 
bination of the four Minkowski functionals. In this sense 
the four Minkowski functionals supply a complete and 
unique characterization of global morphology in three 
dimensions. Table ^displays their direct relation to mor- 
phological quantities known from differential geometry. 
For a smooth body K, they are given as surface inte- 
grals of functions of the principal curvature radii R\, 
R 2 : 



A = 



H 



dA, 

OK 



dA 

dK 



(3) 



X 



^ JdK RlR2 



(for a rigorous derivation of these relations see Mecke, 
Buchert & Wagner 1994). A good overview and exten- 
sive references on integral geometry are given in an ar- 
ticle by Weil (1983) or in the recent book by Schnei- 
der (1993). For further details on Minkowski functionals 
in the cosmological context see Mecke (1994), Buchert 
(1995), Platzoder & Buchert (1995), Schmalzing, Ker- 
scher & Buchert (1995) and especially Mecke, Buchert 
& Wagner (1994). 



4.2 Boundary correction 

For Minkowski functionals there is a concise way of deal- 
ing with boundaries (Mecke & Wagner 1991). Let D be 
the window (the sample geometry) through which we 
look at N clusters. A r — (J. =1 B r (i) is the union of 
balls B r (i) of radius r centered on the i-th cluster, re- 
spectively. In order to obtain a precise quantitative mea- 
sure of the boundary contribution we have to calculate 
the Minkowski functionals M u {A r n D) of the intersec- 
tion of the union of all balls with the window, and the 
Minkowski functionals M M (.D) of the window itself, as 
illustrated in Figure ^. The quantities (A r n D) are 
well suited for the analysis of redshift catalogues but we 
can go further. 

Following Mecke & Wagner (1991) and Schmalzing, 
Kerscher & Buchert (1995), we can extract the volume 
densities of the Minkowski functionals m M (^4 r ) from the 
catalogue if the window contains a fair sample of an er- 
godic and stationary point process. The boundary con- 




Figure 4. Two-dimensional cut through the geometry of 
one part. The shaded area is the set A r D D, in this case D 
is the overall sample geometry (thin line) shrinked by r. 



tribution of the window is then completely removed by 
applying the following recursive formulay 



m M (A r ) = 

M^(A r n D) 
M (D) 



E 

i/=0 



rn v (A r 



M^-AD) (4) 
M (D) ' 



Fava & Santalo (1979) give a mathematically rigorous 
derivation of this formula. In the case of periodic bound- 
ary conditions we have M V {D) = for ^6 {1,2,3} 
and therefore m ti (A r ) = M u (A r n D)/M {D) for 
fi 6 {0, 1, 2, 3}. Generally, for /i = 0, the formula is triv- 
ial. We get for the volume density mo(A r ) (and for the 
void probability function, see below) 



m (A r ) = 



M (A r nD) 
M (D) 



(5) 



Thus, for the removal of boundary contributions of the 
Euler characteristic (and hence the genus, see below) 
in a window D with arbitrary non-periodic boundaries 
one needs to know all the Minkowski functionals. Coles, 
Davies & Pearson (1996) proposed a method for cal- 
culating the genus of isodensity surfaces based on the 
Morse theorem. Their boundary correction is exact for 
periodic cubes. 



4.3 Connection to other statistics 

The most notable Minkowski functional is the Euler 
characteristic \ whose investigation has become a stan- 
dard method of cosmology known as genus-statistics 
(see e.g. Gott, Weinberg & Melott 1987, Melott 1990). 
In three dimensions, for a single body B the Euler char- 
acteristic x(^) is related to the genus g of the surface 
dB by 



X(B) 



\ X (dB) = hl- 



<?)- 



(6) 



2^ ' 2 

The genus is usually calculated from the integral Gaus- 
sian curvature of smoothed isodensity surfaces. The con- 

* we use the convention Y" 1 . x„ = for j < i 
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struction of the isodensity surface involves two para- 
meters, the smoothing length and the density thresh- 
old. In our model we employ the radius as a single diag- 
nostic scale parameter. For simulation data, with peri- 
odic boundaries, the construction of the density field is 
straightforward; for redshift surveys one has to rely on 
boundary corrections (Rhoads, Gott & Postman 1994). 
Since we do not need the density field in our analysis 
we are not concerned with these additional (empirical) 
corrections. 

Another well known statistic is the void probability 
Po (see White 1979) which is directly related to the 
volume density mo according to 



l--Po(r) =m (A r ). 



(7) 



Similar to the connection of the void probability with 
the hierarchy of correlation functions (see Stratonovich 
1963 and White 1979), one finds analogous expres- 
sions for all Minkowski functionals (Mecke 1994, Mecke, 
Buchert & Wagner 1994). 

There is a rule of thumb stating that the first zero 
of the Euler characteristic serves as an estimate for the 
percolation threshold (Mecke & Wagner 1991). 

The dependence of the volume Mo(A r ) on the 
radius r allows the calculation of the Minkowski- 
Bouligand dimension Dm 



lim 



log(M (A r )) 
log(r) 



(8) 



which is equal to the capacity (box counting) dimen- 
sion, giving an upper limit to the Hausdorff dimension 
(see Falconer 1990). As in the case of the ordinary box- 
counting method, any estimate of Dm based on a finite 
number of points is affected by discreteness problems, 
as discussed by Borgani et al. (1993) and Dubrulle & 
Lachieze-Rey (1994). 



4.4 Analysis of the Abell/ACO catalogue 

In order to calculate the densities of Minkowski func- 
tionals from the Abell/ACO sample, we analyzed the 
northern and southern part separately. Then we calcu- 
lated the average between these two parts; since they 
do not intersect, this gives the same results as analyzing 
both parts together. In Figure ^ we sketch the geometry 
of one part. At each radius r, the walls of the original 
window are shrunk by the radius r. All clusters within 
the original window contribute, but clusters lying out- 
side the original window do not. We can then apply the 
boundary removal described in Subsection |l| for each ra- 
dius separately to recover the densities m M (A r ) of the 
Minkowski functionals. Figures ^, ^ and (?) display these 
densities for the various processes described above. 

According to their definition, the densities of the 
Minkowski functionals are expressed in the following 
units: 



[m ] = 1, 

[m 2 ] = (ft _1 Mpc)- 



[mi] = (h^Mpc)" 1 , 
[m 3 ] = (/r'Mpc)- 3 . 



In Figure q we plot the densities of the Minkowski 
functionals for both the Abell/ACO sample (solid lines) 
and a Poisson process (shaded area) . Although the mean 




80 40 
h-'Mpc] 



Figure 5. Densities of the Minkowski functionals for the 
Abell/ACO (average of northern and southern parts) and a 
Poisson process (shaded area) with the same number density. 
The shaded area gives the statistical variance of the Poisson 
process calculated from 100 different realizations. 



values for the Poisson process are known analytically 
(Mecke & Wagner 1991), in Figure |H] we preferred to 
compute numerically the mean values and the standard 
errors of the Minkowski functionals for 100 realizations 
of such a Poisson process within the sample geometry. 

The most prominent feature of all four Minkowski 
functionals are the broader extrema for the Abell/ACO 
data as compared to the results for the Poisson process. 
This is a first indication for enhanced clustering. Let us 
now look at each functional in detail. 

The Minkowski functional mo measures the density 
of the covered volume. On scales between 25ft _1 Mpc 
and 40 ft, - Mpc, mo as a function of r lies slightly below 
the Poisson data. The volume density is lower because of 
the clumping of clusters on those scales. It is remarkable 
that this behaviour occurs in the spatial region where 
clusters are weakly clumped as measured by the con- 
ventional two-point autocorrelation function (the region 
where < £(r) < 1). The Minkowski functionals involve 
correlation functions of every order (see Section |l.3| ) and 
therefore are more sensitive to enhanced clumping. 

The Minkowski functional mi measures the surface 
density of the coverage. It has a maximum at about 
20 h~ Mpc both for the Poisson sample and for the clus- 
ter data. This maximum is due to the granular struc- 
ture of the union set on the relevant scales. At the same 
scales, we find the maximum deviation from the Pois- 
son distribution. The lower values of the cluster data 
mi with respect to the Poisson samples is again the sig- 
nature for the presence of significant clumping of clus- 
ters at these scales. The functional mi shows also a 
positive deviation from the Poisson samples on scales 
(35. . .50) h _1 Mpc where more coherent structures form 
in the union set than in the Poisson samples, keeping 
the surface density larger. 
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The Minkowki functionals 777,2 and 777,3 characterize 
in more detail the kind of spatial coverage provided by 
the union set of balls in the data sample. The density of 
the total mean curvature 7712 of the data reaches a max- 
imum at about 10/i _1 Mpc produced by the dominance 
of convex (positive 777,2) structures. The density 777,2 at 
the maximum is reduced with respect to the Poisson 
sample to about 70% (or more than 3 standard devia- 
tions). The integral mean curvature mi has a zero at a 
scale of 25 /i -1 Mpc (almost the scale of maximum of mi) 
corresponding to the turning-point between structures 
with mainly convex and concave boundaries (negative 
mi). Significant deviations from the Poisson distribu- 
tion occur between this turning point and 40/i _1 Mpc 
due to the smaller mean curvature of the union set of 
the data, probably caused by the interconnection of the 
void regions in the cluster distribution. 

The density of the Euler characteristic 7773 de- 
scribes the global topology of the cluster distribution. 
On small scales all balls are separated. Therefore, each 
ball gives a contribution of unity to the Euler char- 
acteristic and m-}, is proportional to the cluster num- 
ber density. As the radius increases, more and more 
balls overlap and 777,3 decreases. At a scale of about 
20ft _1 Mpc it drops below zero due to the emergence 
of tunnels in the union set (a double torus has x — ~ !)• 
The positive maximum for the Poisson process at scales 
~ 40 h~ Mpc is the signature for the presence of cavi- 
ties. The nearly linear decrease of the Euler characteris- 
tic for the Abell/ACO sample indicates strong cluster- 
ing on scales i; 15 fo -1 Mpc. This confirms the results of 
Bahcall (1988), stating that superclusters consist mainly 
of pairs or triples of strongly correlated (rich) clusters. 
The lack of a significant positive maximum after the 
minimum shows that only a few cavities form. This sug- 
gests on such scales a support dimension for the distri- 
bution of clusters of less than three, (see also Borgani 
et al. 1993); note that the formation of cavities result- 
ing in positive contributions to the Euler characteristic 
is not possible for a support dimension < 2. The pres- 
ence of voids on scales of 30 to 45fo -1 Mpc is shown 
by the enhanced surface area mi and the reduced inte- 
gral mean curvature 7712, while on these scales the Euler 
characteristic 7773 is approximately zero. This does not 
confirm a pure shell model, proposed by Bahcall (1988) 
and used by Mecke, Buchert & Wagner (1994) as the 
"double Poisson process" ; rather, the shape of the clus- 
ter distribution is more likely described by a mixture 
of cavities and tunnels. However, since the scatter in 
the Minkowski functionals on scales above 45 ft -1 Mpc 
is quite high, this last interpretation remains to be con- 
firmed by larger future data sets. 

4.5 Analysis of mock catalogues 

Since the construction of mock samples by the algorithm 
described in Section ^ takes into account all selection ef- 
fects of the Abell/ACO catalogue, we can analyze them 
by using the same procedure as for the real data. The 
considerable deviations of the mock samples' number 
density from the mean cluster number density, as de- 
fined within the whole simulation box, are clearly seen 
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Figure 6. Densities of the Minkowski functionals for the 
Abell/ACO (solid line in both panels) compared to the 
SCDM (shaded area in top panel), and the TCDM (shaded 
area in bottom panel). The shaded area gives la— errorbars 
of the variance among different realizations including cosmic 
variance, as explained in the text. 



in the scatter of the density of the Euler characteris- 
tic 777,3 at the smallest radii, which is proportional to 
the number density (see Figures ^ and ^ . This suggests 
that even a sample as large as the Abell/ACO cata- 
logue is not a 'fair sample' as discussed by Buchert and 
Martinez (1993), since it is still significantly influenced 
by cosmic variance. 

A comparison of the four SCDM (or the three 
ACDM) realizations shows that this scatter in the num- 
ber density is the main source for the variance seen in 
the Minkowski functionals. The statistical variance in- 
troduced by the scatter of different realizations of the 
initial density field is always negligible if compared to 
the observer-to-observer scatter in each box, associated 
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Figure 7. Densities of the Minkowski functionals for the 
Abell/ACO (solid line in both panels) compared to the 
ACDM (shaded area in top panel), and the BSI (shaded area 
in bottom panel). The shaded area gives la— errorbars of the 
variance among different realizations including cosmic vari- 
ance, as explained in the text. 



with the number density fluctuations between different 
mock samples. In this sense the error bands shown in 
the plots include cosmic variance; only the 'robustness' 
(a result of the additivity property) and the high signif- 
icance of the Minkowski functionals still allow for dis- 
crimination between the different Dark Matter cosmolo- 
gies. 

In Figure ^ we compare the densities of the func- 
tionals for SCDM and TCDM models to the Abell/ACO 
ones. Both models show too little clustering on small 
scales, as it is clearly seen by the enhanced maxima of 
the surface area mi and the integral mean curvature mi, 
as well as in the flatter decrease of the Euler character- 
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Figure 8. Differences of the densities of Minkowski func- 
tionals for the SCDM (solid line), TCDM (dotted), ACDM 
(long dashed), and the BSI (short dashed). 



istic vri3. Additionally, the higher volume mo indicates 
weak clumping also on large scales. 

Figure fj] displays the same quantities as Figure |^ 
but for the ACDM and BSI models. Although small de- 
viations from the Abell/ACO results are still present, 
again pointing towards weak clumping, both models 
perform much better than SCDM and TCDM. 

This is also confirmed by the plot in Figure ^, where 
we show the differences between the functionals of the 
Abell/ACO data, m^ ata , and the functionals of the mock 
samples, m™° ck . Again we see that ACDM and BSI de- 
scribe the cluster distribution much better than SCDM 
and TCDM. At scales above 45/i _1 Mpc the differences 
between the Minkowski functionals of the Abell/ACO 
and the mock samples are dominated by random fluc- 
tuations in the Abell/ACO. 

In order to quantify the statistical significance of 
these deviations we calculate the risks R^, with an er- 
ror weighted quadratic cost function, for all functionals 
of all models versus those of the Abell/ACO sample. 
This risk is equal to the x 2 distance used in maximum 
likelihood analysis (see e.g. Frieden 1991). For N r radii 
it is given by 



c (r0 



(9) 



where <7™° (r) is the standard error of the Minkowski 
functionals at radius r over several mock catalogues, 



\r) = v/(m™ ock (r) 2 ) - (r 



<(r)>* 



(10) 



In Figure ^| we show the i? M values for the four con- 
sidered DM models. After excluding r = and r > 
45 h~ Mpc, we limit the analysis to = 18 radial bins. 
Therefore, a constant departure of la should result in 
a risk — 18, which is shown for reference in Fig- 
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Figure 9. The risk values for simulations versus Abell/ACO 
for the volume measure i?o (empty boxes), for the surface 
density Ri (crosses), for the mean curvature R2 (stars), and 
for the Euler characteristic i?3 (filled boxes). 



rather than isolated voids. Our results from the mock 
catalogues agree with the expectation that the SCDM 
and TCDM models do not describe the degree of clus- 
ter clumping as inferred from the data in the scale range 
(15. . . 45) h~ Mpc. Moreover the formation of cavity ag- 
gregates is not reproduced convincingly. On the other 
hand, the ACDM and the BSI models provide a reason- 
ably good although not perfect description of the data. 

The main source of uncertainty in the analysis 
arises from the fluctuations in the number density of 
objects among the mock catalogues. The deviations of 
the model results from the data are on the la level. 
Nevertheless we think that the results are reliable, since 
they are inferred from the behaviour of our morpholog- 
ical measures over an extended range of the diagnostic 
length scale parameter r. 

The family of Minkowski functionals focuses on 
global features of spatial patterns and includes, besides 
a topological descriptor for the connectivity, also ge- 
ometrical measures for the size and shape. Thus the 
present approach provides a unifying frame for the 
analysis of cosmic structures which comprises the void 
probability function as well as the genus statistics; our 
method is introduced not to replace but to complement 
traditional tools such as low-order correlation functions. 



ure |9j as a dotted line. SCDM shows a significant depar- 
ture from Abell/ACO, as a result of the reduced clus- 
tering. The differences of TCDM are even larger than 
those of SCDM (cf. also Figure ^). The higher risk val- 
ues for TCDM may be caused by underestimating 
(we recall that we analyzed only one TCDM simulation, 
i.e. eight mock catalogues, but four SCDM simulation, 
i.e. 32 mock catalogues). As for ACDM and BSI, they 
are confirmed to fit the data much better; the difference 
with respect to Abell/ACO data is below the lcr level 
for all the functionals, but showing also a tendency to- 
wards too weak clumping. The volume density (i.e. the 
void probability) alone is not sufficiently discriminating 
(see also Figures ^, , while the integral mean curvature 
rri2 is even more selective than the Euler characteristic 
m 3 . 



The numerical code used for calculating the Minkowski 
functionals is available and can be obtained by sending 
e-mail to buchert@stat.physik.uni-muenchen.de. 
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5 SUMMARY AND CONCLUSION 

We have calculated scale-dependent Minkowski func- 
tionals for spatial patterns induced by the point set of an 
Abell/ACO cluster redshift sample and compared these 
functionals with those obtained from a number of mock 
cluster catalogues obtained from different cosmological 
models. We find significant deviations in the morpholog- 
ical features of the observed cluster clumping from the 
Poisson distribution on scales (15. . . 50) /i _1 Mpc, which 
have not been explored previously in such detail. In par- 
ticular, in the range (15. . . 45) /i _1 Mpc the results indi- 
cate a clumping of clusters on a support with dimension 
less than three (see also Borgani et al. 1994). On scales 
(25. . .45) h~ 1 Mpc the behavior of the Minkowski func- 
tionals mi (area), 7772 (integral mean curvature) and 7773 
(Euler characteristic) may be interpreted as a cluster ag- 
gregation exhibiting cavities and interconnected tunnels 
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